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Abstract 

We investigate the possibility to discriminate between different pairs of CP 
non-conserving mechanisms inducing the neutrinoless double beta (/3/3)oi/-decay by 
using data on (/3/3)o,,-decay half-lives of nuclei with largely different nuclear ma- 
trix elements (NMEs). The mechanisms studied are: light Majorana neutrino ex- 
change, heavy left-handed (LH) and heavy right-handed (RH) Majorana neutrino 
exchanges, lepton charge non-conserving couplings in SUSY theories with imparity 
breaking giving rise to the "dominant gluino exchange" and the "squark-neutrino" 
mechanisms. The nuclei considered are 76 Ge, 82 Se, 100 Mo, 130 Te and 136 Xe. Four 
sets of nuclear matrix elements (NMEs) of the decays of these five nuclei, derived 
within the Self-consistent Renormalized Quasiparticle Random Phase Approxima- 
tion (SRQRPA), were employed in our analysis. While for each of the five single 
mechanisms discussed, the NMEs for 76 Ge, 82 Se, 100 Mo and 130 Te differ relatively 
little, the relative difference between the NMEs of any two nuclei not exceeding 
10%, the NMEs for 136 Xe differ significantly from those of 76 Ge, 82 5e, 100 Mo and 
130 Te, being by a factor ~ (1.3 — 2.5) smaller. This allows, in principle, to draw 
conclusions about the pair of non-interfering (interfering) mechanisms possibly in- 
ducing the (/3/3)o;/-decay from data on the half- lives of 136 Xe and of at least one 
(two) more isotope(s) which can be, e.g., any of the four, 76 Ge, 82 Se, 100 Mo and 
130 Te. Depending on the sets of mechanisms considered, the conclusion can be in- 
dependent of, or can depend on, the NMEs used in the analysis. The implications 
of the EXO lower bound on the half-life of 136 Xe for the problem studied are also 
exploited. 

1 Introduction 

If neutrinoless double beta ( / 9 / 9)o I /-decay will be observed, it will be of fundamental im- 
portance to determine the mechanism which induces the decay. In [T] we have considered 

1 Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 
Sofia, Bulgaria 
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the possibility of several different mechanisms contributing to the neutrinoless double 
beta (/3/3)ot/-decay amplitude in the general case of CP nonconservation [§. The mecha- 
nisms discussed are the "standard" light Majorana neutrino exchange, exchange of heavy 
Majorana neutrinos coupled to (V-A) currents, exchange of heavy right-handed (RH) 
Majorana neutrinos coupled to (V+A) currents, lepton charge non- conserving couplings 
in SUSY theories with i?-parity breaking H. Of the latter we have concentrated on the 
so-called "dominant gluino exchange" and "squark-neutrino" mechanisms. Each of these 
mechanisms is characterized by a specific fundamental lepton number violating (LNV) 
parameter r] K , the index k labeling the mechanism. The parameter r] K will be complex, 
in general, if the mechanism k does not conserve the CP symmetry. In [HI UJ [8] the 
authors analised the possibility to identify the mechanisms of (/3/3)o^-decay using data 
on the half-lives of several isotopes. The indicated (and other) specific mechanisms were 
considered, assuming that only one of these mechanisms is triggering the decay @. In [I] 
we have investigated in detail the cases of two "non-interfering" [5] and two "interfering" 
mechanisms generating the ( J 5 ( 5)o 1 /-decay H, using as input hypothetical (/3/3)oj/-decay half- 
lives of the three isotopes 76 Ge, 100 Mo and 130 Te. Four sets of nuclear matrix elements 
(NMEs) of the decays of these three nuclei, derived within the Self-consistent Renor- 
malized Quasiparticle Random Phase Approximation (SRQRPA) [TTj 112], were utilized: 
they were calculated in [T] with two different nucleon-nucleon (NN) potentials (CD-Bonn 
and Argonne), "large" size single-particle spaces and for two values of the axial coupling 
constant qa = 1.25; 1.0. 

If the (/3/3)oy-decay is induced by two non-interfering mechanisms, which for con- 
creteness were considered in [T] to be [5] the light left-handed (LH) and the heavy RH 
Majorana neutrino exchanges, one can determine the squares of the absolute values of 
the two LNV parameters, characterizing these mechanisms, |?7„| 2 and |?7r| 2 , from data 
on the half-lives of two nuclear isotopes. This was done in [T| using as input all three 
possible pairs of half-lives of 76 Ge, 100 Mo and 130 Te, chosen from specific intervals and 
satisfying the existing experimental constraints. It was found that if the half-life of one 
of the three nuclei is measured, the requirement that \rj u \ 2 > and |?7_r| 2 > ("posi- 
tivity condition") constrains the other two half-lives (and the ( / 9 / 9) 0! y-decay half-life of 
any other (/3/3)o!/-decaying isotope for that matter) to lie in specific intervals, determined 
by the measured half-life and the relevant NMEs and phase-space factors. This feature 
is common to all cases of two non-interfering mechanisms generating the (/3/5)o^-decay 

2 The case of two CP conserving mechanisms generating the (/3/3)oy-decay was considered in [5]. 

3 For a more detailed description of these mechanisms and references to the articles where they were 
originally proposed see pQ and, e.g., OH]. 

4 In \W\ the authors derived a general effective (/3/3)oi,-decay Lagrangian without specifying the 
mechanisms generating the different terms in the Lagrangian. Assuming that only one of the terms in 
the Lagrangian is operative in (/3/3)o l /-decay and using the existing lower limit on the half-life of 76 Ge, 
in |10| constraints on the effective CP conserving couplings constant multiplying the different terms in 
the Lagrangian were obtained. 

5 Two mechanisms contributing to the (/3/3)o^-decay amplitude can be non-interfering or interfering, 
depending on whether their interference term present in the (/3/3)o y -decay rate is suppressed (and negligi- 
ble) or not. In the case of the five mechanisms considered in [T] and listed above, all pairs of mechanisms, 
which include the exchange of the heavy RH Majorana neutrinos coupled to (V+A) currents as one of 
the mechanisms, can be shown to be of the non-interfering type, while those that do not involve the 
exchange of heavy RH Majorana neutrinos belong to the interfering class. 
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The indicated specific half-life intervals for the various isotopes were shown to be stable 
with respect to the change of the NMEs (within the sets of NMEs employed) used to 
derive them. The intervals depend, in general, on the type of the two non-interfering 
mechanisms assumed to cause the (/?/3)oi/-decay. However, these differences in the cases 
of all possible pairs of non-interfering mechanisms considered were found to be extremely 
small. Using the indicated difference to get information about the specific pair of non- 
interfering mechanisms possibly operative in (/3/3) 0j ,-decay requires, in the cases studied in 
[T|, an extremely high precision in the measurement of the (/3/3) j,-decay half-lives of the 
isotopes considered ( 76 Ge, 100 Mo and 130 Te), as well as an exceedingly small uncertainties 
in the knowledge of the (/3/3)o^-decay NMEs of these isotopes. The levels of precision 
required seem impossible to achieve in the foreseeable future. One of the consequences 
of this results is that if it will be possible to rule out one pair of the considered in [1] 
non-interfering mechanisms as the cause of (/3/9)o^-decay, most likely one will be able to 
rule out all of them. 

The dependence of the physical solutions for 1 77^ | 2 and \t]r\ 2 obtained on the NMEs 
used, was also studied in [T|. It was found that the solutions can exhibit a significant, 
or a relatively small, variation with the NMEs employed, depending on the hypothetical 
values of the half-lives of the two isotopes utilized as input for obtaining the solutions. 
This conclusion is valid for all other pairs of non-interfering mechanisms considered in 
[T]. In the case when two interfering mechanisms are responsible for the (/3/3)o^-decay, 
the squares of the absolute values of the two relevant parameters and the interference 
term parameter, which involves the cosine of an unknown relative phase a of the two 
fundamental parameters, can be uniquely determined, in principle, from data on the 
half-lives of three nuclei. We have analyzed in [T] in detail the case of light Majorana 
neutrino exchange and gluino exchange. In this case the parameters which are determined 
from data on the half-lives are \r] u \ 2 , |?7a'| 2 , Vx being that of the gluino exchange, and 
z = 2 cos a \r) u \ \r)y\. The physical solutions for these parameters have to satisfy the 
conditions \r] u \ 2 > 0, |?7a'| 2 > and — 2|r/ 1/ ||r/y | < z < 2 ] 77^ | ] | . The latter condition 
implies that given the half-lives of two isotopes, T\ and T 2 , the half- life of any third 
isotope T 3 is constrained to lie is a specific interval, if the mechanisms considered are 
indeed generating the (/3/3)o;,-decay. If further the half-life of one isotope T\ is known, 
for the interference to be constructive (destructive), the half-lives of any other pair of 
isotopes T 2 and T 3 , should belong to specific intervals. These intervals depend on whether 
the interference between the two contributions in the (/3/3)oi,-decay rate is constructive 
or destructive. We have derived in [I] in analytic form the general conditions for i) 
constructive interference (z > 0), ii) destructive interference (z < 0), iii) |r/^| 2 = 0, 
|?7a'| 2 7^ ; iv) \t] u \ 2 7^ 0, |?7a'| 2 = and v) z — 0, \r] u \ 2 ^ 0, \r]y\ 2 ^ 0. We have found 
that, given T 1; a constructive interference is possible only if T 2 lies in a relatively narrow 
interval and T 3 has a value in extremely narrow intervals. Numerically the intervals for 
T 2 and T 3 are very similar to the intervals one obtains in the case of two non-interfering 
mechanisms (within the set considered in [T]). The intervals of values of T 2 and T 3 
corresponding to destructive interference are very different from those corresponding 
to the cases of constructive interference and of the two non-interfering (/3/3)oi/-decay 
mechanisms we have considered. Within the set of (/3/3) 0j ,-decay mechanisms studied 
by us, this difference can allow to discriminate experimentally between the possibilities 
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of the (/3/3)o l /-decay being triggered by two " destructively interfering" mechanisms or by 
two "constructively interfering" or by two non-interfering mechanisms. 

The "degeneracy" of the predictions of the pairs of non-interfering mechanisms of 
(/3/3)oi/-decay for the interval of values of the half-life of a second nucleus, given the half-life 
of a different one from the three considered, 76 Ge, 100 Mo and 130 Te, is a direct consequence 
of a specific property of the NMEs of the three nuclei. Namely, for each of the five 
single mechanisms discussed (the light LH, heavy LH and heavy RH Majorana neutrino 
exchanges, the gluino exchange and the "squark-neutrino" mechanism), the NMEs for the 
three nuclei differ relatively little, the relative difference between the NMEs of any two 
nuclei not exceeding 10% B \M i)K - M jtK \/(0.5\M ijK + M jjK \) £ 0.1, i ^ j = 1,2,3 = 76 
Ge, 100 Mo, 130 Te. This feature of the NMEs of 76 Ge,' 100 Mo and 130 Te makes it impossible 
to discriminate experimentally not only between the four different pairs of non-interfering 
mechanisms, considered in [T], using data on the half- lives of 76 Ge, 100 Mo and 130 Te, 
but also between any of these four pairs and pairs of interfering mechanisms when the 
interference is constructive. The indicated "degeneracy" of the predictions of different 
pairs of mechanisms possibly active in (/3/3)oi,-decay can be lifted if one uses as input 
the half-lives of nuclei having largely different NMEs. One example of such a nucleus is 
136 Xe, whose NMEs for the five mechanisms studied in [T], as we are going to show, differ 
significantly from those of 76 Ge, 82 Se, 100 Mo and 130 Te. 

In the present article we investigate the potential of combining data on the half-lives 
of im Xe and of one or more of the four nuclei 76 Ge, 82 Se, 100 Mo and 130 Te, for discriminat- 
ing between different pairs of non-interfering or interfering mechanisms of (/3/3)oi,-decay. 
We consider the same five basic mechanisms used in the study performed in [T], namely, 
the "standard" light Majorana neutrino exchange, exchange of heavy Majorana neutrinos 
coupled to (V-A) currents, exchange of heavy right-handed (RH) Majorana neutrinos cou- 
pled to (V+A) currents, dominant gluino exchange and the squark-neutrino mechanism. 
The last two are related to lepton charge non- conserving couplings in SUSY theories with 
i?-parity breaking. 

The paper is organized as follows. In Section 2 we give a brief overview on the possible 
mechanisms that can induce the (/3/3) 0! ,-decay considered in this work. In Section 3 we 
analyze the case of two mechanisms active in the decay in the non-interfering and in the 
interfering regime, combining the recent experimental results reported by EXO, while the 
Section 4 contains our conclusions. 

2 Mechanisms of (/3/3)oz/-Decay Considered and Nuclear 
Matrix Elements Employed 

The mechanisms of (/3/3)o;/-decay we are going to consider in the present article are the 
same five mechanisms considered in pQ: i) the light Majorana neutrino exchange, char- 
acterized by the dimensionless LNV parameter r] u = <m>/m e , <m> and and m e being 
the ( / 9 / 9)o I /-decay effective Majorana mass (see, e.g., 011]) and the electron mass, respec- 
tively; ii) exchange of heavy "left-handed" (LH) Majorana neutrinos coupled to (V-A) 

6 The general implications of this "degeneracy" of the NMEs of 76 Ge, 100 Mo and 130 Te for testing the 
mechanisms under discussion in the case of CP invariance were investigated in [T3] . 
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currents, characterized by the LNV parameter tjl) in) exchange of heavy "right-handed" 
(RH) Majorana neutrinos coupled to (V+A) currents, characterized by the LNV param- 
eter rjR. We consider also two possible mechanisms in SUSY theories with R-parity 
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Figure 1: Feynman diagrams for the (/3/3)o;,-decay, generated by the light and heavy LH 
Majorana neutrino exchange (left panel) and the heavy RH Majorana neutrino exchange 
(right panel). 
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Figure 2: Feynman diagrams for (/3/3) 0[ /-decay due to the gluino exchange mechanism. 

non- conservation: iv) the "gluino exchange" mechanisms, characterized by the LNV pa- 
rameter rjx', and v) the "squark-neutrino" mechanism, characterized by the parameter 
ri u _ q . All five LNV parameters we have introduced above r] v , T]l, Vr> Vx' an< ^ Vu-q are 
effective. The dependence of a given LNV effective parameter on the fundamental param- 
eters of the theory in which the corresponding mechanism is possible/arises, is discussed 
m detail in [I] and we are not going to repeat this discussion here. Instead we will 
give just a graphical representation of each of the five mechanisms in terms of the corre- 
sponding leading order Feynman diagrams. For the light Majorana neutrino, heavy LH 
Majorana neutrino (N)~l) and heavy RH Majorana neutrino (A^) exchange mechanisms 
they are given in Fig. [TJ for the gluino exchange and squark-neutrino mechanisms we 
show them in Figs. [2] and [3j 

7 An extensive list of references to the original articles is also given in pQ . 
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Figure 3: Feynman diagrams for (/3/3) 0! ,-decay due to the squark-neutrino mechanism at 
the quark-level |14j . 

In the present analysis we employ four sets of nuclear matrix elements (NMEs) of the 
decays of the five nuclei of interest, 76 Ge, 82 Se, 100 Mo, 130 Te and 136 Xe, derived within the 
Self-consistent Renormalized Quasiparticle Random Phase Approximation (SRQRPA) 
[TT j [T2" ] . The SRQRPA takes into account the Pauli exclusion principle and conserves the 
mean particle number in correlated ground state. 

For each of the five nuclei, two choices of single-particle basis are considered. The 
intermediate size model space has 12 levels (oscillator shells N=2-4) for 76 Ge and 82 Se, 
16 levels (oscillator shells N=2-4 plus the f+h orbits from N=5) for 100 Mo and 18 levels 
(oscillator shells N=3,4 plus f+h+p orbits from N=5) for 130 Te and 136 Xe. The large 
size single particle space contains 21 levels (oscillator shells N=0-5) for 76 Ge, 82 Se and 
100 Mo, and 23 levels for 130 Te and 136 Xe (N=l-5 and i orbits from N=6). In comparison 
with previous studies |15], we omitted the small space model which is not sufficient to 
describe realistically the tensor part of the (/S/^o^-decay nuclear matrix elements. 

The single particle energies were obtained by using a Coulomb-corrected Woods- 
Saxon potential. Two-body G-matrix elements we derived from the Argonne and the 
Charge Dependent Bonn (CD-Bonn) one-boson exchange potential within the Brueckner 
theory. The schematic pairing interactions have been adjusted to fit the empirical pairing 
gaps |16j . The particle-particle and particle- hole channels of the G-matrix interaction of 
the nuclear Hamiltonian H are renormalized by introducing the parameters g pp and g p h, 
respectively. The calculations have been carried out for g p h = 1.0. The particle-particle 
strength parameter g pp of the SRQRPA is fixed by the data on the two-neutrino double 
beta decays [151 HZ]- in the calculation of the (/3/3)oiy-decay NMEs, the two-nucleon short- 
range correlations derived from same potential as residual interactions, namely from the 
Argonne or CD-Bonn potentials, were considered |18| . 

The calculated NMEs W °", M'^, M'°y and M'f are listed in Table [U For 76 Ge, 
82 Se, 100 Mo and 130 Te they are taken from ref. [T|, while for 136 Xe the results are new 
|19] . We note that these NMEs are significantly smaller (by a factor 1.3 - 2.5) when 
compared with those for 76 Ge, 82 Se and 100 Mo. The reduction of the Oz//3/3-decay NMEs 
of the 136 Xe is explained by the closed neutron shell for this nucleus. A sharper Fermi 
surface leads to a reduction of this transition. This effect is clearly seen also in the case 
of M J of double magic nucleus 48 Ca |19] . 
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By glancing the Table [T] we see that a significant source of uncertainty is the value of 
the axial-vector coupling constant and especially in the case of matrix elements M'y 
and M'^ u . Further, the NMEs associated with heavy neutrino exchange are sensitive also 
to the choice of the NN interaction, the CD-Bonn or Argonne potential. These types of 
realistic NN interaction differ mostly by the description of the short-range interactions. 
Although in Table [T] we present results for NMEs of the nuclei of interest, calculated 
using both medium and large size single particle spaces within the SRQRPA method, in 
the numerical examples we are going to present further we will use the NMEs for a given 
nucleus, with the large size single particle space in both cases of Argonne and CD-Bonn 
potentials and for = 1.25; 1.00 (i.e., altogether four NMEs). 

Table 1: Nuclear matrix elements M'°J (light neutrino mass mechanism), M' ^ (heavy 
neutrino mass mechanism), M'°y (trilinear R-parity breaking SUSY mechanism) and 
M'°q (squark mixing mechanism) for the Ozv/3/3-decays of 76 Ge, 100 Se, 100 Mo, 130 Te and 
136 Xe within the Selfconsistent Renormalized Quasiparticle Random Phase Approxima- 
tion (SRQRPA). G 0u (E , Z) is the phase-space factor. We notice that all NMEs given in 
Table [T]are real and positive. The nuclear radius is R = 1.1 fm A 1 / 3 . 
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3 (/?/3)oi/- Decay Induced by Two Mechanisms 

The observation of (j3 (3) 0u -dec&y of several different isotopes is crucial for obtaining in- 
formation about the mechanism or mechanisms that induce the decay. In the analysis we 
are going to perform we will employ the lower bound obtained by the EXO collaboration 
on the (/3/3) 0l ,-decay half-life of 136 Xe [20]: 

T^ 2 { im Xe) > 1.6 x 10 25 y (90%CL). (1) 

We use also the lower limits on the (/3/3)o;,-decay half-lives of Ge, 82 Se and 100 Mo, and 
of 130 Te reported by the Heidelberg-Moscow [23], NEM03 [22] and CUORICINO [23] 
experiments, respectively, as well as as well the 76 Ge half-life reported in [21] (see also 
J): 

Ty 2 C 6 Ge) > 1.9 x 10 25 y [21], Ty 2 ( 82 Se) > 3.6 x 10 23 y 



T^i^Mo) > 1.1 x 10 24 y [22], T^^Te) > 3.0 x 10 24 y [23] . (2) 

•76 

/2\ 



Z&^Ge) = 2.23±°; 44 x 10 25 y 



Following [T], we will consider two cases: 

1. ( ( 5/3)o I y-decay induced by two mechanisms whose interference term in the (/3{3)o v - 
decay rate is negligible |f| ; 

2. (/3/3)oi/-decay triggered by two CP non- conserving mechanisms whose interference 
term cannot be neglected. 

In the case 1, given the two mechanisms A and B, the inverse of the (/3/3)ojy-decay 
half-life for a given isotope (Aj, Zj) reads: 

\Va\ 2 \M%\ 2 + \ VB \*\M%\ 2 , (3) 



TiGi 



where the index i denotes the isotope. The values of the phase space factor G® U (E, Z), 
and of the NMEs M'^ A and M'® U B for the mechanisms we will consider and for the isotopes 
76 Ge, S2 Se, 100 Mo, 130 Te and 136 Xe of interest, are listed in Table[U The LNV parameters 
are defined in section [21 If the two mechanisms A and B inducing the decay are interfering 
and CP non-conserving (case 2), the inverse of the (/3/3)oj,-decay half-life of the isotope 
(Ai, Zi) can be written as: 

|^| 2 |Mt A | 2 + |^| 2 |Mt s | 2 + 2cosa|Mt A ||MtBll^lkB|. (4) 



T^G^(E,Z) 



Here a is the relative phase of t]a and t]b- 

If the (/3/3)oi/-decay is caused by two non-interfering mechanisms, the LNV parameters 
\tia\ 2 and \tib\ 2 characterizing the mechanisms, can be determined, in principle, from 



8 This possibility is realized when, e.g., the electron currents (responsible for the emission of the two 
electrons in the final state), associated with the two mechanisms considered, have opposite chiralities. 
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data on the (/3/3)o^-decay half-lives of two isotopes, i.e., by solving a system of two linear 
equations. In the case of two interfering CP non- conserving mechanisms, the values of 
the two parameters |?7a| 2 and \tIb\ 2 and of the cosine of the relative phase a of 7] a and 
T]b can be obtained from data on the (/3/3)o ! y-decay half-lives of three isotopes, i.e., by 
solving a system of three linear equations. 

As was noticed and discussed in detail in [I], a very important role in identifying 
the physical solutions for |?7a| 2 an d \^b\ 2 of the corresponding systems of two or three 
equations is played by the "positivity conditions" |?7a| 2 > and \t]b\ 2 > 0. If one of 
the two mechanisms inducing the (/3/3)o;,-decay is the "standard" light neutrino exchange, 
additional important constraint on the positive solutions of the relevant systems of two or 
three linear equations can be provided by the upper limit on the absolute neutrino mass 
scale set by the Moscow and Mainz 3 if /3-decay experiments [26| [27] : m{v e ) < 2.3 eV. 
In the case of (/3/3) 0; ,-decay, this limit implies a similar limit on the effective Majorana 
masa^l |<m>| < 2.3 eV. The latter inequality translates into the following upper bound 
on the corresponding LNV dimensionless parameter \rj u \ 2 = (|<m>| /m e ) 2 : 

1^1 2 x 10 12 < 21.2. (5) 

The KATRIN 3 if /3-decay experiment |27j . which is under preparation, is planned to have 
approximately a 10 times better sensitivity to m(u e ) than that achieved in the Moscow 
and Mainz experiments. If the designed sensitivity limit of |<m>| < 0.2 eV (90% CL.) 
will be obtained in the KATRIN experiment, it would imply the following rather stringent 
upper limit on \r] v \ 2 : 

\r] v \ 2 x 10 12 < 0.16. (6) 

In this work we will derive numerical results using the NMEs calculated with the large 
size single particle basis ("large basis") and the Argonne potential ("Argonne NMEs"). 
We report also results obtained with NMEs calculated with the Charge Dependent Bonn 
(CD-Bonn) potential ("CD-Bonn NMEs") and compared them with those derived with 
the Argonne NMEs. 



3.1 Two Non-interfering Mechanisms 

In this case the solutions for the corresponding two LNV parameters \t]a\ 2 and |?7_b| 2 
obtained from data on the (/3/3)o„-decay half-lives of the two isotopes (A iy Z^) and (A,-, Zj), 
are given by [1]: 

\VA\ = i TTTjZ ioi , Tl n v , i » r,()v W\ TJrtD n > \Vb\ 







WW, - 




\ 2 /TiGi 


\Mf A 


21 




2 ~\M%\ 


2 \M'Xa\ 2 



(7) 

It follows from eq. (jTJ that [T] if one of the two half-lives, say T i; is fixed, the positivity 
conditions |/?a| 2 > and \t]b\ 2 > can be satisfied only if Tj lies in a specific "positivity 



We recall that for mi. 2, 3 > 0.1 eV the neutrino mass spectrum is quasi-degenerate (QD), mi = rri2 
2 > AmlijIAml! 



m.3 = m, uij 3> Ato^ij |Ato3x|. In this case we have m(v e ) = m and |<m>| < m. 
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interval". Choosing for convenience always Aj < A{ we get for the positivity interval [T]: 



d \M% 2 d Mf A 2 

Gj \M%\ 2 l ~ ~ Gj \M'f A \ 2 l ' y ' 

where we have used | Mf A | 2 / 1 M'f A | 2 > \M%\ 2 / \M'^ B \ 2 . In the case of \M%\ 2 /\M'% A \ 2 < 
\M'°i B \ 2 /\M' 2 h ' B \ 2 , the interval of values of Tj under discussion is given by: 



r I M' 0u I 2 r I M' 0u I 2 

Gj. \ m i,A\ rp < T < ^£ i,Bl T /q\ 
Gj |M'°^| 2 *" j ~ G j \Mf B \i l ' 



Condition ([8]) is fulfilled, for instance, if A is the heavy right-handed (RH) Majorana 
neutrino exchange and B is the light Majorana neutrino exchange in the case of Argonne 
NMEs (see Table [I]). The inequality in eq. ([8]) (or 0) has to be combined with the 
experimental lower bounds on the half-lives of the considered nuclei, T^ in . If, e.g., 
^Tmin * s the lower bound of interest for the isotope (Ai,Zi), i.e., if T > T^ in , we get 
from eq. ([8]): 

r \M' 0u I 2 

rp > I l,B\ rpexp / 1r) N 

~~ Gj \M'°j" B \ 2 min ' 1 j 

The lower limit in eq. (TTU]) can be larger than the existing experimental lower bound on 
Tj. Indeed, suppose that T = T^ 2 ( U6 Xe), Tj = T^ 2 ( 76 Ge) and that the (/3/3) 0l/ -decay 
is due by the standard light neutrino exchange and the heavy RH Majorana neutrino 
exchange. In this case the positivity conditions for \rj u \ 2 and |?7r| 2 imply for the Argonne 
and CD-Bonn NMEs corresponding to g A = 1.25 (1.0): 

T 0v ( 7e Ge) 

1.90 (1.85) < T y 2 ( i3 6Xe) < 2-70 (2.64), (Argonne NMEs) ; (11) 
T 0u ( 76 Ge) 

1-30 (1.16) < T ^ 2 (136Xe) < 2.47 (2.30), (CD-Bonn NMEs) . (12) 
Using the EXO result, eq. (pQ), and the Argonne NMEs we get the lower bound on 

1/2 ^ 



Tx/ 2 C 6 Ge) > 3.03 (2.95) x 10 25 y. (13) 

This lower bound is significantly bigger that the experimental lower bound on Ty 2 ( 7G Ge) 
quoted in eq. (T5]). If we use instead the CD-Bonn NMEs, the limit we obtain is close to 
the experimental lower bound on T^ 2 ( 7e Ge): 

Tx/ 2 {™Ge) > 2.08 (1.85) x 10 25 y. (14) 

For illustrative purposes we show in Fig. H] the solutions of equation ([Tj) for \r) u \ 2 
and \t] R \ 2 derived by fixing T^ 2 ( 76 Ge) to the best fit value claimed in [23, T^ 2 ( 76 Ge) = 
2.23 x 10 25 (see eq. ([2])). As Fig. H] shows, the positive (physical) solutions obtained 
using the Argonne NMEs are incompatible with the EXO result, eq. ([T|), and under the 
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The values of \rj u \ 2 (solid lines) and |?7r| 2 (dashed lines) obtained for 



as a function of Ty 2 ( 136 Xe), using the Argonne (left 



Figure 4: 

T°J 2 ( 76 Ge) = 2.23 x 10^ y g 
panel) and CD-Bonn (right panel) NMEs corresponding to qa = 1-25 (thick lines) and 
gA = 1 (thin lines). The region of physical (positive) solutions for qa = 1-25 are delimited 
by the two vertical lines. The solid horizontal line corresponds to the prospective upper 
limit from the KATRIN experiment [27], while the thick solid vertical line indicates the 
EXO lower bound [2D]. The gray areas correspond to excluded values of \r] v \ 2 and 1 77^ | 2 . 



assumptions made and according to our oversimplified analysis, are ruled out. At the 
same time, the physical solutions obtained using the CD-Bonn NMEs are compatible 
with the EXO limit for values of \r] u \ 2 and \t)r\ 2 lying in a relatively narrow interval. 

We consider next a second example of two non-interfering (/3/3) ^-decay mechanisms, 
i.e., (/9/3)o^-decay induced by the heavy RH Majorana neutrino exchange and the gluino 
exchange. Setting, as above, T« = TfY 2 ( 13S Xe) and Tj = Ty^^Ge), we get for the posi- 
tivity intervals using the Argonne or CD-Bonn NMEs corresponding to ^ = 1.25 (1.0): 

T 0v ( 76 Ge) 

2-70 (2.64) < T y 2 (136Xe) < 2-78 (2.67), (Argonne NMEs) ; (15) 
T 0v ( 76 Ge) 

1-30 (1.16) < r ^ 2 (136Xe) < 4.43 (4.25), (CD-Bonn NMEs) , (16) 

The lower bound on T®J 2 ( 7e Ge) following from the EXO limit in the case of the Argonne 
NMEs obtained with g A = 1.25 (1.0) reads: 

Tff 2 ( n Ge) > 4.31 (4.22) x 10 25 y. (17) 

This lower limit is by a factor of 2.27 bigger than the experimental lower limit quoted in 
eq. (j2J). It is also incompatible with the 4a range of values of Ty 2 ( 76 Ge) found in [24] . 
The lower bound obtained using the CD-Bonn NMEs is less stringent: 

Tx/ 2 {™Ge) > 2.08 (1.85) x 10 25 y. (18) 

In Fig. |we show that if T^ 2 ( 7e Ge) = 2.23 x 10 25 , the recent EXO lower limit allows 
positive (physical) solutions for the corresponding two LNV parameters |?7a'| 2 and |?7r| 2 
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only for the CD-Bonn NMEs and for values of |?7a'| 2 and |?7_r| 2 lying in a very narrow 
interval. 



76 Ge (r 1/2 =2.23 10 25 ) and 136 Xe 



76 Ge (r 1/2 =2.23 10 25 ) and 136 Xe 




d_ 10 



7.5 X 10 24 



8.X10 24 8.5 Xl0 2 " 

T V i( 136 Xe) [y] 



9. X 10 2 




5.X10 24 



1.x 10 25 1.5 xlO 25 

ri /2 ( 136 Xe) [y] 



2. X 10 2S 



Figure 5: The same as in Fig. H] but for the values of the rescaled parameters |?7a'| 2 
(solid lines) and |?7r| 2 (dashed lines). 

We get similar results for the third pair of non-interfering mechanisms - the squark- 
neutrino exchange and the heavy RH neutrino exchange. Indeed, using the Argonne 
NMEs corresponding to qa = 1-25 (1.0) we find for the positivity interval: 



2.52 (2.40) < 



T^Xe) 



< 2.70 (2.64) , 



The EXO lower bound in this case implies: 



Ti/ 2 (™Ge) > 4.03 (3.84) x 10 25 y. 
From the NMEs computed with the CD-Bonn potential we get 



1.30(1.16) < 



T°; 2 (i36 Xe ) 



< 2.95 (2.8i; 



and 



TxU™Ge) > 2.08 (1.85) x 10 25 y- 



(19) 



(20) 



(21) 



(22) 



In the case the non-interfering LH and RH heavy Majorana neutrino exchanges, the 
NMEs for the two mechanisms coincide and the system of equation in (EJ) reduces to a 
relation between the half-lives of the two considered isotopes: 



T T G,\M%? 
3 l G 3 \M'%? 



(23) 



In this case the EXO lower bound implies the following lower limits on Tj = Ty 2 ( 76 Ge) 
for the sets of NMEs we are considering for g& = 1.25 (1.0): 

Ty 2 ( 76 Ge) > 4.31 (4.22) x 10 25 y (Argonne NMEs) , (24) 
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T^ 2 ( 7G Ge) > 2.08 (1.85) x 10 25 y (CD-Bonn NMEs). (25) 

The range of positive solutions for the LNV parameters in equation flH]) shifts towards 
larger values if Tj is increased. As we noticed in PQ, if the experimentally determined 
interval of allowed values of the ratio Tj/Ti of the half-lives of the two isotopes considered, 
including all relevant uncertainties, lies outside the range of positive solutions for |//a| 2 
and \tib\ 2 i one would be led to conclude that the (/3/3)o I /-decay is not generated by the 
two mechanisms under discussion. 

Assuming the half-lives of two isotopes, say, of 76 Ge and 136 Xe, Ty 2 ( w Ge) = T\ and 
Ty 2 ( Xe) = T 2 , to be known, and (/3/3)o;,-decay triggered by a pair of non-interfering 
mechanisms A and B, one can always use the physical solutions for \t)^b \ 2 (Ti, ^2), 
obtained using the two half-lives Ti i2 (in eq. ([7])), to find the range of the half-life of a 
third isotope: 

^ = G 3 (|^(T 1 ,T 2 )| 2 |M'^r + |^(T 1 ,T 2 )| 2 |M'^r), (26) 

In Tables [2] and |3] we give numerical predictions based on this observation. Fixing the 
half-life of 76 Ge to T\ = 10 26 y and assuming the l36 Xe half-life T 2 lies in an interval 
compatible with the existing constraints, the system of two equations is solved and the 
values of \t]a\ 2 > and \tib\ 2 > thus obtained are used to get predictions for the half-life 
of a third isotope, in this case 82 Se, 100 Mo and 130 Te. The mechanisms considered are 
a) light and heavy RH Majorana neutrino exchanges (Table |2]) and b) gluino and heavy 
RH Majorana neutrino exchanges (Table [3]). It follows from the results shown in Tables 
[2] and |3] that the intervals of allowed values of the half-lives of 82 Se, 100 Mo and 130 Te thus 
obtained i) are rather narrow 0, and ii) exhibit weak dependence on the NMEs used to 
derive them (within the sets of NMEs considered) . 

One can use eq. (126]) and the lower bound, e.g., on T 1 ° / i/ 2 ( 136 Xe) reported by the EXO 
experiment, to derive a lower bound on one of the half-lives involved in the study of two 
non- interfering mechanisms, say T\. Indeed, we can set T3 = Ty 2 ( 136 Xe), in eq. (1261) . use 
the explicit form of the solutions for \t)a F b v \ 2 (Ti, T 2 ) and apply the existing EXO lower 
bound. We get: 

1 Di D 2 1 

_ — — J < 27) 

T 3 NTx NT 2 1.6xl0 25 y' V ; 

where 

Dl = Ii {\ M 'V\M'tB? ~ \M't A \ 2 \M%\ 2 ) , D 2 = g [\M'f A \ 2 \M%\ 2 - \M'X A \ 2 \M%\ 2 
and 

N = \M'° 2 U A \ 2 \M%\ 2 - \M'\ V A \ 2 \M%\ 2 . (28) 
Using further the positivity constraint given in eq. ([S]), 

aTi < T 2 < bTi , (29) 



10 We note that the experimental lower bounds quoted in eq. ([2]) have to be taken into account since, 
in principle, they can further constrain the range of allowed values of 1 77^4 1 2 and \t]b\ 2 and of the half-life 
of the third isotope of interest. 
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Table 2: Predictions using Argonne and CD-Bonn NMEs corresponding to ^=1.25 (5U=1 
in parenthesis) in the case of two non-interfering mechanism: light and heavy RH Ma- 
jorana neutrino exchanges. The physical solutions for \rj u \ 2 and |?7_r| 2 derived for given 
half-lives of 76 Ge and 136 Xe, are used to obtain predictions for the half-lives of 82 Se, 
W0 Mo and 130 Te. The 76 Ge half-life was set to T( 76 Ge) = 10 26 yr, while the interval of 
values of the 13S Xe half-life was determine from the positivity conditions. 



Argonne NMEs 




Positive solutions 


Predictions 




3.71(3.79) ■ 


2.30(2.34) 

10 25 < T( 136 Xe) < 5.27(5.42) • 10 25 1.45(1.46) ■ 

1.76(1.78) ■ 


• 10 25 < T( S2 Se) < 2.39(2.49) ■ 
10 25 < T( 100 Mo) < 1.80(1.76) 
10 25 < T( 13Q Te) < 2.44(2.49) 


10 25 
• 10 25 
10 25 


CD-Bonn NMEs 




Positive solutions 


Predictions 




4.04(4.35) ■ 


2.30(2.33) 

10 25 < T( 136 Xe) < 7.71(8.63) • 10 25 1.44(1.45) • 

1.65(1.68) ■ 


• 10 25 < T( S2 Se) < 2.39(2.48) ■ 
10 25 < T( 100 Mo) < 1.78(1.74) 
10 25 < T( 13Q Te) < 2.21(2.27) 


10 25 

.1 25 

10 25 



Table 3: The same as in Table [2] but for the gluino and RH heavy Majorana neutrino 
exchange mechanisms. 



Argonne NMEs 




Positive solutions 


Predictions 


3.60(3.74) ■ 


2.27(2.32) 

10 25 < T( 136 Xe) < 3.71(3.79) • 10 25 1.43(1.459) • 

1.76(1.78) ■ 


• 10 25 < T( S2 Se) < 2.30(2.34) • 10 25 
10 25 < T( 100 A/o) < 1.45(1.460) • 10 25 
10 25 < T( 130 Te) < 1.80(1.86) • 10 25 


CD-Bonn NMEs 




Positive solutions 


Predictions 


2.26(2.35) ■ 


2.27(2.32) 

10 25 < T( 136 Xe) < 7.71(8.63) • 10 25 1.43(1.4542) • 

1.65(1.68) ■ 


• 10 25 < T( S2 Se) < 2.30(2.33) • 10 25 
10 25 < T( lm Mo) < 1.44(1.4543) • 10 25 
10 25 < T( 130 Te) < 1.75(1.81) • 10 25 
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where b = \M'°{ A \ 2 /\M'° 2 " A \ 2 > a = \M%\ 2 / \M%\ 2 , we get the following lower limit 
on T\ from eq. (I2"7j) : 

T,>T, >1.6xl0 25 y (30) 

1 - 3 \N bNJ y \N bNJ K ' 

This lower bound on 2\ depends only on T 3 and on the NMEs |M'^| 2 and \M'°j U B \ 2 , 
i,j = 1,2,3. We give examples of predictions based on the eq. (I3"0"j) in Tables H] and 
We notice that, for the NMEs used in the present study, the EXO lower limit on 
Ty 2 ( l36 Xe) sets a lower bound on the half-lives of the other isotopes considered by us 
that usually exceed their respective current experimental lower bounds. 

Table 4: Lower bound on Ti from eq. (I3"0"j) using the EXO limit on Ty 2 ( 136 Xe), eq. (CQ), 
and the Argonne and CD-Bonn NMEs corresponding to ^=1.25 (^=1), in the case of 
two non-interfering mechanisms - light and heavy RH Majorana neutrino exchanges. See 
text for details. 



m,r 2 ) 


Argonne g^=1.25 (1.0) 


CD-Bonn g A =1.25 (1.0) 


130 Te - 76 Ge 
100 Mo - 76 Ge 
82 Se - 76 Ge 


T( 130 Te)> 7.40 (7.35)- 10 24 
T( 100 Mo)> 5.45 (5.19)- 10 24 
T( 82 Se)> 7.25 (7.37)- 10 24 


T( 130 Te)> 3.43 (3.11) ■ 10 24 
T( 130 Te)> 3.00 (2.70) • 10 24 
T( 82 Se)> 4.76 (4.32) • 10 24 



Table 5: The same as in Table H] for the gluino and heavy RH Majorana neutrino 
exchange mechanisms. See text for details. 



m,T 2 ) 


Argonne g^=1.25 (1.0) 


CD-Bonn g A =l-25 (1.0) 


130 Te - 76 Ge 
100 Mo - 76 Ge 
82 Se - 76 Ge 


T( 130 Te)> 7.59 (7.53)- 10 24 
T( 100 Mo)> 6.25 (6.16)- 10 24 
T( 82 Se)> 9.90 (9.87)- 10 24 


T( 13U Te)> 3.43 (3.11) • 10 24 
T( 130 Te)> 3.00 (2.70) • 10 24 
T( 82 Se)> 4.76 (4.32) • 10 24 



3.2 Discriminating between Different Pairs of Non-interfering 
Mechanisms 

The first thing to notice is that, as it follows from Table 1, for each of the four dif- 
ferent mechanisms of (/3/3)o^-decay considered, the relative difference between NMEs 
of the decays of w Ge, 82 Se, 100 Mo and 130 Te does not exceed approximately 10%: 
(M'f x - Mf x )/(0.5(M'^ X + Mf x )) & 0.1, where i^j = 76 Ge, 82 5e, 100 Mo, 130 Te, and 
X denotes any one of the four mechanisms discussed. As was shown in [I], this leads 
to degeneracies between the positivity intervals of values of the ratio of the half-lives 
of any two given of the indicated four isotopes, corresponding to the different pairs of 
mechanisms inducing the (/3/3)o^-decay The degeneracies in question make it practically 

11 The inequality in eq. was derived assuming that D^/N > 0. In the case of D2/N < one has 
to interchange a and b in it. 
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impossible to distinguish between the different pairs of (/3/3)o v -decay mechanisms, con- 
sidered in pP| and in the present article, using data on the half-lives of two or more of the 
four nuclei 76 Ge, 82 Se, 100 Mo and 130 Te. At the same time, it is possible, in principle, to 
exclude them all using data on the half-lives of at least two of the indicated four nuclei 

IB- 

In contrast, the NMEs for the (PP)o u -decaj of 136 Xe, corresponding to each of the 
four different mechanisms we are considering are by a factor of ~ (1.3 — 2.5) smaller than 
the (/3/3)o„-decay NMEs of the other four isotopes listed above: {M'f x - M'° iX ) / M'^ x ) = 
(0.3 - 1.5), where i = 136 Xe and j = 76 Ge, 82 Se, 100 Mo, 130 Te (see Figs. Eland [7]). As a 
consequence, using data on the half-life of im Xe as input in determining the positiv- 
ity interval of values of the half-life of any second isotope lifts to a certain degree the 
degeneracy of the positivity intervals corresponding to different pairs of non-interfering 
mechanisms. This allows, in principle, to draw conclusions about the pair of mechanisms 
possibly inducing the (/3/3)o I y-decay from data on the half-lives of 136 Ae and a second 
isotope which can be, e.g., any of the four considered above, 76 Ge, 82 Se, 100 Mo and 
130 Te. 

To be more specific, it follows from eqs. f JTT]) . f fT5|) . f fl9|) . (1231) and Table [1] that if the 
Argonne (CD-Bonn) NMEs derived for = 1-25 (1.0) are correct, all four pairs of mech- 
anisms of (/3/3)o!/-decay discussed by us will be disfavored, or ruled out, if it is established 
experimentally that R( 76 Ge, 136 Xe) > 2.8 (4.5) or that R( 76 Ge, 136 Xe) < 1.8 (1.1), where 
R( 76 Ge, 136 Xe) = T°Y 2 ( 76 Ge)/T°Y 2 ( 136 Xe). Further, assuming the validity of the Argonne 
NMEs, one would conclude that the light and heavy RH Majorana neutrino exchanges 
are the only possible pair of mechanisms operative in (P(3) 0u -decay if it is found experi- 
mentally that 1.9 < R( 76 Ge, 13G Xe) < 2.4. For 1.9 < R( 76 Ge, 136 Xe) < 2.6, i) the gluino 
and RH Majorana neutrino exchanges, and ii) the LH and RH heavy Majorana neutrino 
exchanges, will be disfavored or ruled out. One finds similar results using the CD-Bonn 
NMEs. The numbers we quote in this paragraph should be considered as illustrative only. 
In a realistic analysis one has to take into account the various relevant experimental and 
theoretical uncertainties. 

We analyze next the possibility to discriminate between two pairs of non-interfering 
mechanisms triggering the (/3/3)o!/-decay when the pairs share one mechanism. Given 
three different non-interfering mechanisms A, B and C, we can test the hypothesis of 
the (/3/3)o;,-decay induced by the pairs i) A + B or ii) C + B, using the half-lives of the 
same two isotopes. As a consequence of the fact that B is common to both pairs of 
mechanisms, the numerators of the expressions for |?7a| 2 and 1 770 1 2 , as it follows from eq. 
([7]), coincide. Correspondingly, using the half-lives of the same two isotopes would allow 
us to distinguish, in principle, between the cases i) and ii) if the denominators in the 
expressions for the solutions for |?7a| 2 and \rjc\ 2 have opposite signs. Indeed, in this case 
the physical solutions for \tia\ 2 in the case i) and \r)c\ 2 in the case ii) will lie either in the 
positivity intervals ([HD and respectively, or in the intervals (JH]) and (JS]). Thus, the 
positivity solution intervals for |?7a| 2 and \r]c\ 2 would not overlap, except for the point 
corresponding to a value of the second isotope half-life where r/A = rjc = 0. This would 
allow, in principle, to discriminate between the two considered pairs of mechanisms. 

It follows from the preceding discussion that in order to be possible to discriminate 
between the pairs A + B and C + B of non-interfering mechanisms of (/3/3)oi/-decay, the 



16 



following condition has to be fulfilled: 



Det | 


(\Mf A \ 2 


\M%\ 2 ) 




2 \M%\ 2 - 


\M% 


2 \M'f A \ 2 


Det | 


r\M%\ 2 

K\M'tc\ 2 


\M%?\ 
\M%?) 


\M% 


2 \M%\ 2 - 


\M% 


2 \M%\ 2 



<0. (31) 



This condition is satisfied if one of the following two sets of inequalities holds: 

M% - M% M% ~ M% M% - M% 

n M% < M% < m ,o :a , w 

rn M% - Mf A M% - M% M% - M% 
] Mf A < M% < M% ' 1 ^ 

One example of a possible application of the preceding results is provided by the 
mechanisms of light Majorana neutrino exchange (A), RH heavy Majorana neutrino 
exchange (B) and gluino exchange (C) and the Argonne NMEs. We are interested in 
studying cases involving 136 Xe since, as it was already discussed earlier, the NMEs of 
136 Xe differ significantly from those of the lighter isotopes such as 7Q Ge (see Tabled]). 
Indeed, as can be shown, it is possible, in principle, to discriminate between the two pairs 
A + B and C + B of the three mechanisms indicated above if we combine data on the 
half-life of 136 Xe with those on the half-life of one of the four isotopes 76 Ge, 82 Se, 100 Mo 
and 130 Te, and use the Argonne NMEs in the analysis. In this case the inequalities (13"2|) 
are realized, as can be seen in Fig. |6l where we plot the relative differences (M'^ v — 
M'") I 'M' " for the Argonne NMEs where the indices % and j refer respectively to 136 Xe 
and to one of the four isotopes 76 Ge, 82 Se, 100 Mo and 130 Te. In the case of the CD-Bonn 
NMEs (Fig. [7]), the inequalities fl32l) or fl33|) do not hold for the pairs of mechanisms 
considered. The inequalities given in eq. (132|) hold, as it follows from Fig. [TJ if, e.g., the 
mechanisms A, B and C are respectively the heavy RH Majorana neutrino exchange, the 
light Majorana neutrino exchange and the gluino exchange. 

The preceding considerations are illustrated graphically in Figs. |S] and In Fig. 
[8]we use Ti = T^ 2 ( 76 Ge) and Tj = Tty 2 ( 136 Xe) and the Argonne (left panel) and CD- 
Bonn (right panel) NMEs for the decays of 76 Ge and 136 Xe to show the possibility of 
discriminating between the two pairs of non-interfering mechanisms considered earlier: 
i) light Majorana neutrino exchange and heavy RH Majorana neutrino exchange (RHN) 
and ii) heavy RH Majorana neutrino exchange and gluino exchange. The 76 Ge half-life 
is set to Ti = 5 x 10 25 y, while that of 136 Xe, Tj, is allowed to vary in a certain interval. 
The solutions for the three LNV parameters corresponding to the three mechanisms 
considered, \r] u \ 2 , |?7r| 2 and |?7a'| 2 > obtained for the chosen value of Tj and interval of 
values of Tj, are shown as functions of Tj. As is clearly seen in the left panel of Fig. [HJ if 
1 77^ 1 2 , |?7r| 2 and |?7a'| 2 are obtained using the Argonne NMEs, the intervals of values of Tj 
for which one obtains the physical positive solutions for \r] u \ 2 and |??a'| 2 ) do not overlap. 
This makes it possible, in principle, to determine which of the two pairs of mechanisms 
considered (if any) is inducing the (/3/3) 0; y-decay The same result does not hold if one 
uses the CD-Bonn NMEs in the analysis, as is illustrated in the right panel of Fig. [HJ In 
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NMEs relative difference using Xe and £4=1.25 



NMEs relative difference using 136 Xe and gA=l 




76Q,. 82 S( , 100 M() 130 Tf . 



Figure 6: The relative differences between the Argonne NMEs (MJ - M'i)/M'i, 
where i= im Xe and j = 76 Ge, 82 Se, 100 Mo, 130 Te, for g A = 1.25 (left panel) and g A = 1 
(right panel) and for three different non-interfering mechanisms: light Majorana neutrino 
exchange (circles), RH heavy Majorana neutrino exchange (squares) and gluino exchange 
(diamonds). See text for details. 



NMEs relative difference using 136 Xe and g A =1.25 NMEs relative difference using 136 Xe and gA=l 




' 6 Ge 82 Se ™Mo 130 Te 



"Ge 82 Se "Mo I30 Te 



Figure 7: The same as in Fig. E]for the CD-Bonn NMEs. See text for details. 



18 



this case none of the inequalities (132]) and ( 133]) is fulfilled, the intervals of values of Tj 
for which one obtains physical solutions for \r] u \ 2 and \f]\\ 2 overlap and the discrimination 
between the two pairs of mechanisms is problematic. 

We show in Fig. |9] that the features of the solutions for \r] v \ 2 and |?7a| 2 we have 
discussed above, which are related to the values of the relevant NMEs, do not change if 
one uses in the analysis the half-lives and NMEs of 136 Xe and of another lighter isotope 
instead of 76 Ge, namely, of 100 Mo. 



Argonne Potential Cd-Bonn Potential 




r 1/2 ( 136 Xe) [y] T 1/2 ( 136 Xe) [y] 

Figure 8: Solutions for the LNV parameters corresponding to two pairs of non-interfering 
mechanisms: i) \ri v \ 2 and \^r\ 2 (dot-dashed and dashed lines) and ii) |?7a'| 2 and \ f]R\ 2 (solid 
and dotted lines). The solutions are obtained by fixing Tj = T®Y 2 (J 6 Ge) = 5 x 10 25 y and 
letting free Tj = Ty 2 ( 136 Xe) and using the sets of Argonne (left panel) and CD-Bonn 
(right panel) NMEs calculated for g& = 1.25 (thick lines) and g& = 1 (thin lines). The 
range of positive solutions in the case of Argonne NMEs and g& = 1.25 is delimited by 
the two vertical dashed lines. The horizontal solid line corresponds to the prospective 
upper limit |<m>| < 0.2 eV |27| . The thick solid vertical line indicates the EXO lower 
limit on T°Y 2 ( 136 Xe) [20] . Sec text for details. 



3.3 Two Interfering Mechanisms 

We analyze in the present Section the possibility of „-decay induced by two inter- 
fering CP-non-conserving mechanisms. This case is characterized by three parameters: 
the absolute values and the relative phase of the two LNV parameters associated with 
the two mechanisms. They can be determined, in principle, from data on the half-lives of 
three isotopes, T i: % — 1, 2, 3. Given 7\ 2,3 and denoting by A and B the two mechanisms, 
one can set a system of three linear equations in three unknowns, the solution of which 
reads [T|: 

\Va\ 2 = -^, \Vb\ 2 = -^-, z = 2 cos a\r] A \\r] B \ = —£ , (34) 
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Figure 9: Solutions for the LNV parameters of two pairs of non-interfering (/3/3) „-decay 
mechanisms i) |?7„| 2 and |?7_r| 2 (dot-dashed and dashed lines) and ii) |t7a'| 2 an d \vr\ 2 (solid 
and dotted lines) obtained by fixing = T 1 ° / y 2 ( 100 Mo) = 6.5 x 10 24 yr and letting free 
Tj = T°Y 2 ( 136 Xe). The other notations are the same as in Fig. |8j See text for details. 



where D, Di, Dj and D k are the following determinants: 



D 



{M'X A y 



{M%f 



B) 



M%M% 
Mf B M'f A 
Mt B M'f A 



i/TiPi 
i/T k G k 



Mf B Mf A 
M'f B M'f A 
Mf B Mf A 



Di 



i/TiGi 
i/TjGj 
i/T k G k 



{M%f 



b) 



(MTa) 2 (M'Z 



M%M'X A 
Mf B Mf A 



BJ 



i/TjG, 
t/T k G k 



(35) 



(36) 

As in the case of two non-interfering mechanisms, the LNV parameters must be non- 
negative | r] A | 2 > 00 and \tib\ 2 > 0, and in addition the interference term must satisfy the 
following condition: 



- 2\r) A \\r] B \ < 2 cos a\f] A \\ri B \ < 2\r] A \\r) B \ 



(37) 



These conditions will be called from here on "positivity conditions". 

Using the positivity conditions it is possible to determine the interval of positive 
solutions for one of the three half-life, e.g., T k , if the values of the other two half-lives 
in the equations have been measured and are known. The condition on the interference 
term in equation ([8]) can considerably reduce the interval of values of T k where \f) A \ 2 > 
and \f] B \ 2 > 0. In Table [6] we give examples of the constraints on T k following from the 
positivity conditions for three different pairs of interfering mechanisms: light Majorana 
neutrino and supersymmetric gluino exchange; light Majorana neutrino exchange and 
heavy LH Majorana neutrino exchange; gluino exchange and heavy LH Majorana neutrino 
exchange. It follows from the results shown in Table El in particular, that when T( 76 Ge) 
is set to T( 76 Ge) = 2.23 x 10 25 ; 10 26 y, but T( 130 Te) is close to the current experimental 
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lower limit, the positivity constraint intervals of values of T( 136 Xe) for the each of the 
three pairs of interfering mechanisms considered are incompatible with the EXO lower 
bound on T( 136 Xe), eq. jQ. 



Table 6: Ranges of the half-live of 136 Xe for different fixed values of the half-lives of 76 Ge 
and 130 Te in the case of three pairs of interfering mechanisms: light Majorana neutrino 
exchange and gluino exchange (upper table); light Majorana and heavy LH Majorana 
neutrino exchanges (middle table); gluino exchange and and heavy LH Majorana neutrino 
exchange (lower table). The results shown are obtained with the "large basis" = 1.25 
Argonne NMEs. One star (two stars) indicate that the EXO bound constrains further 
(rules out) the corresponding solution. 



T?y y ] (fixed) 


T?; 2 [y] (fixed) 


Allowed Range 


T(Ge) = 2.23 ■ 10 25 ** 
T(Ge) = 10 26 ** 
T(Ge) = 2.23 • 10 25 
T(Ge) = 10 26 


T(Te) = 3 • 10 24 
T(Te) = 3 • 10 24 
T(Te) = 3 • 10 25 
T(Te) = 3 • 10 25 


2.95 • 10 24 < T(Xe) < 5.65 • 10 24 
3.43 • 10 24 < T(Xe) < 4.66 ■ 10 24 
1.74 • 10 25 < T(Xe) < 1.66 • 10 26 
2.58 • 10 25 < T(Xe) < 6.90 • 10 25 



T?; 2 [y] (fixed) 


T?; 2 [y] (fixed) 


Allowed Range 


T(Ge) = 2.23 • 10 25 ** 
T(Ge) = 10 26 ** 
T(Ge) = 2.23 ■ 10 25 
T(Ge) = 10 26 


T(Te) = 3 • 10 24 
T(Te) = 3 • 10 24 
T(Te) = 3 ■ 10 25 
T(Te) = 3 ■ 10 25 


4.93 • 10 24 < T(Xe) < 6.21 • 10 24 
5.23 • 10 24 < T(Xe) < 5.83 ■ 10 24 
3.95 • 10 25 < T(Xe) < 8.25 ■ 10 25 
4.68 • 10 25 < T(Xe) < 6.61 ■ 10 25 



T?; 2 [y] (fixed) 


T^ 2 [y] (fixed) 


Allowed Range 


T(Ge) = 2.23 • 10 25 ** 
T(Ge) = 10 26 * 
T(Ge) = 2.23 • 10 25 ** 
T(Ge) = 10 26 * 


T(Te) = 3 • 10 24 
T(Te) = 3 • 10 24 
T(Te) = 3 • 10 25 
T(Te) = 3 • 10 25 


5.59 ■ 10 23 < T(Xe) < 1.26 • 10 25 
1.21 ■ 10 24 < T(Xe) < 4.71 • 10 25 
1.05 ■ 10 24 < T(Xe) < 2.42 • 10 24 
3.32 • 10 24 < T(Xe) < 2.16 • 10 25 



We consider next a case in which the half-life of 136 Xe is one of the two half-lives 
assumed to have been experimentally determined. The (/3/3)o-decay is supposed to be 
triggered by light Majorana neutrino and gluino exchange mechanisms with LFV param- 
eters \r] u \ 2 and |?7a'| 2 - We use in the analysis the half-lives of 
which will be denoted for simplicity respectively as T 1; T 2 and T 3 
bounds on T i: i = 1,2,3, given in eq. (j2J), are taken into account, the conditions for 



76 Ge, 136 Xe and 130 Te, 
Once the experimental 



destructive interference, i.e., for cos a < 0, are given by: 



z < : < 



1.9 x 10 25 < Ti < 1.90T 2 , T 3 > 
1.90T 2 < T x < 2.78T 2 , T 3 > 
T x > 2.78T 2 , T 3 > 



9.64T!T 2 



16.32Ti + 8.59T 2 ' 
3.82T!T 2 

6.33T! + 3.66T 2 ' 
7.33T!T 2 

11.94Ti + 7.6 1T 2 
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where we have used the "large basis" g& = 1.25 Argonne NMEs (see Table [T]). The 
conditions for constructive interference read: 

9.64T\T 2 3.82T!T 2 
1.90T 2 < Ti < 2.29T 2 , — - < T 3 < 



,vn. 16.32T\ + 8.59T 2 ~ ~ 6.3321 + 3.66T 2 ,o ql 

* > u • S 7.33TiT 2 3.82T!T 2 ^) 

2.29T 2 < Ti < 2.78T 2 , — < T 3 < — . 

{ ' 11.94Ti + 7.61T 2 - 3 - 6.3321 + 3.66T 2 

If we set, e.g., the 76 Ge half-life to the value claimed in |2"1] Ti — 2.23 x 10 25 y, we 
find that only destructive interference between the contributions of the two mechanisms 
considered in the ( / 9/3)o I /-decay rate, is possible. Numerically we get in this case 

3.44T 2 

T-? > . (40) 

3 5.82 + 1.37 x 10- 25 T 2 1 ; 

For 1.37 x 1(T 25 T 2 < 5.82 one finds: 

T( 130 Te) > 0.59 T( 136 Xe) > 9.46 x 10 24 y , (41) 

where the last inequality has been obtained using the EXO lower bound on T( 136 Xe). 
Constructive interference is possible for the pair of interfering mechanisms under discus- 
sion only if T( 76 Ge) > 3.033 x 10 25 y. 

The possibilities of destructive and constructive interference are illustrated in Figs. [TD] 
and [TH respectively. In these figures the physical allowed regions, determined through the 
positivity conditions, correspond to the areas within the two vertical lines (the solutions 
must be compatible also with the existing lower limits given in eq (j2J)). For instance, 
using the Argonne "large basis" NMEs corresponding to = 1-25 and setting T( 76 Ge) = 
2.23 x 10 25 y and T( 130 Te) = 10 25 y, positive solutions are allowed only in the interval 
1.60 x 10 25 < T( 136 Xe) < 2.66 x 10 25 y (Fig. HJJJ. As can be seen in Figs. M and EJ 
a constructive interference is possible only if T 2 = T( 13G Xe) lies in a relatively narrow 
interval and T3 = T( 130 Te) is determined through the conditions in eq. (1391) . 

Next, we would like to illustrate the possibility to distinguish between two pairs of 
interfering mechanisms i) A+B and ii) B+C, which share one mechanism, namely B, 
from the data on the half-lives of three isotopes. In this case we can set two systems 
of three equations, each one in three unknowns. We will denote the corresponding LNV 
parameters as i) |?7a| 2 , \Vb\ 2 and ii) |?7s| 2 and \r]c\ 2 , while the interference parameters 
will be denoted as i) z and ii) z' . Fixing two of the three half-lives, say Tj and Tj, the 
possibility to discriminate between the mechanisms A and C relies on the dependence of 
\t)a\ 2 and \r)c\ 2 on the third half-life, T^. Given T« and Tj, it will be possible to discriminate 
between the mechanisms A and C if the two intervals of values of where | t)a\ 2 > 
and \t]c\ 2 > 0, do not overlap. If, instead, the two intervals partially overlap, complete 
discrimination would be impossible, but there would be a large interval of values of (or 
equivalently, positive solutions values of the LNV parameters) that can be excluded using 
present or future experimental data. In order to have non-overlapping positive solution 
intervals of Tk, corresponding to |?7a| 2 > and \t]c\ 2 > 0, the following inequality must 
hold: 

{M'l A M% ~ M'XaM'Ib^M'IaM'Xb - M%MXb) . . 

{M%M% ~ M'XbM'Xc) (M'XbM'Xc - M%M'Xc) ' 1 j 
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Figure 10: Left panel: the values of \r] u \ 2 x 10 10 (thick solid line) and |?7a'| 2 x 10 14 
(dotted line), obtained as solutions of the system of equations (J3J) for fixed values of 
T( 76 Ge) = 2.23 x 10 25 y and T( 130 Te) = 10 25 y, and letting T^ 2 ( im Xe) free. The physical 
allowed regions correspond to the areas within the two vertical lines. Right panel: the 
values of the phase a in the allowed interval of values of Ty 2 ( 136 Xe), corresponding to 
physical solutions for \rj v \ 2 and \r]\>\ 2 . In this case cos a < and the interference is 
destructive. See text for details. 




1.6 xlO 25 1.65 xlO 25 1.7 xlO 25 1.75 xlO 25 1.6 xlO 25 1.65 xlO 25 Uxltf 

T 1/2 ( 1M Xe) [y] J 1/2 ( 136 Xe)[y] 



Figure 11: Left panel: the same as in Fig. [10] but for T( 76 Ge) = 3.5 x 10 25 y and 
T( 130 Te) = 8.0 x 10 25 y. The interval of values oiTty 2 ( 130 Xe) between i) the vertical solid 
and right dashed lines ii) the two vertical dashed lines, and iii) the vertical solid and left 
dashed lines, correspond respectively to i) physical (non-negative) solutions for \rj„\ 2 and 
\r)y\ 2 , ii) constructive interference (z > 0), and iii) destructive interference (z < 0). Right 
panel: the corresponding values of the phase a as a function of T®" 2 ( 130 Xe). Constructive 
interference is possible only for values of T^ 2 ( l30 Xe) between the two vertical dashed lines. 
See text for details. 



The above condition can be satisfied only for certain sets of isotopes. Obviously, whether 
it is fulfilled or not depends on the values of the relevant NMEs. We will illustrate this on 
the example of an oversimplified analysis involving the light Majorana neutrino exchange, 
the heavy LH Majorana neutrino exchange and the gluino exchange as mechanisms A, 
B and C, respectively, and the half-lives of 76 Ge, 130 Te and 136 le: T x = T( 76 Ge), 
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T 2 = T( 130 Te) and T 3 = T( 136 Xe). Fixing T x = 2.23 x 10 25 y and T 3 = 1.6 x 10 25 y 
(the EXO 90% C.L. lower limit), we obtain the results shown in Fig. H2J As it follows 
from Fig. [121 i n the case of the Argonne NMEs (left panel), it is possible to discriminate 
between the standard light neutrino exchange and the gluino exchange mechanisms: the 
intervals of values of T 2 , where the positive solutions for the LNV parameters of the two 
pairs of interfering mechanisms considered occur, do not overlap. Further, the physical 
solutions for the two LNV parameters related to the gluino mechanism are excluded by 
the CUORICINO limit on T( 130 Te) |23] . This result does not change with the increasing 
of T 3 . Thus, we are lead to conclude that for T 3 > 1.6 x 10 25 y and 7\ given by 
the value claimed in [23], of the two considered pairs of possible interfering (/3/3)o^-decay 
mechanisms, only the light and heavy LH Majorana neutrino exchanges can be generating 
the decay. The solution for \r] u \ 2 must be compatible with the upper limit |<m>| < 2.3 
eV f2E \ |2"7 ] , indicated with a solid horizontal line in Fig. [T2J In the right panel of Fig. [T2] 
we plot also the solutions obtained with the CD-Bonn NMEs. In this case is not possible 
to discriminate between the two considered pair of mechanisms since the condition in eq. 
(142]) is not satisfied. 

Another interesting example is the case in which A is the light Majorana neutrino 
exchange, B is the gluino exchange and C the heavy LH Majorana neutrino exchange, i.e., 
we try to discriminate between i) the light neutrino plus gluino exchange mechanisms, 
and ii) the heavy LH Majorana neutrino plus gluino exchange mechanisms. We fix, like 
in the previous case, the values for T\ = 2.23 x 10 25 y and T 3 = 1.6 x 10 25 y. The results 
of this analysis are plotted in Fig. [13j Since the condition in eq. (142|) is now satisfied for 
NMEs obtained either with the Argonne potential or with the CD-Bonn potential, in this 
case it is possible, in principle, to discriminate between the the two pair of mechanisms 
independently of the set of NMEs used (within the sets considered by us). This result 
does not change with the increasing of T 3 . Hence, as far as 7\ is fixed to the value claimed 
in [21] and the limits in eq. ([2]) are satisfied, the two intervals of values of T 2 , in which 
the "positivity conditions" for i) \rj v \ 2 , \i]y\ 2 and z, and for ii) |?7a'| 2 , IvnI 2 an d z', are 
satisfied, are not overlapping (Fig. [T3"|) . 



4 Conclusions and Summary 

We have investigated the possibility to discriminate between different pairs of CP non- 
conserving mechanisms inducing the neutrinoless double beta (j3f3) o u -decay by using 
data on (/3/3)o^-decay half-lives of nuclei with largely different nuclear matrix elements 
(NMEs). The mechanisms studied are: light Majorana neutrino exchange, heavy left- 
handed (LH) and heavy right-handed (RH) Majorana neutrino exchanges, lepton charge 
non- conserving couplings in SUSY theories with i?-parity breaking giving rise to the 
"dominant gluino exchange" and the "squark-neutrino" mechanisms. Each of these mech- 
anisms is characterized by a specific lepton number violating (LNV) parameter r] K , where 
the index k labels the mechanism. For the five mechanisms listed above we use the nota- 
tions k = v, L, R, A', q, respectively. The parameter rj K will be complex, in general, if the 
mechanism k does not conserve the CP symmetry. The nuclei considered are 76 Ge, 82 Se, 
100 Mo, 130 Te and 136 Xe. Four sets of nuclear matrix elements (NMEs) of the (/3/3) ^-decays 
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Figure 12: The parameters \r] u \ 2 x 10 10 (solid line) and |?7l| 2 X 10 14 (dotted line) of the 
light and heavy LH Majorana neutrino exchange mechanisms, and |?7y| 2 x 10 14 (dashed- 
dotted line) and |?7l| 2 x 10 14 (dashed line) of the gluino and heavy LH Majorana neutrino 
exchange mechanisms, obtained from eq. (j34p using the Argonne NMEs (left panel) and 
CD-Bonn NMEs (right panel), corresponding to Qa = 1-25 (thick lines) and g& — 1 (thin 
lines), for T^ 2 ( w Ge) = 2.23 x 10 25 y T^ 2 ( 136 Xe) = 1.60 x 10 25 y and letting Ty 2 ( 130 Te) 
free. See text for details. 




r 1/2 ( 130 Te) [y] r 1/2 ( 130 Te) [y] 

Figure 13: The same as in Fig. [121 but for i) \r] u \ 2 x 10 10 (thick solid line) and \r)y\ 2 x 
10 14 (thick dotted line) of the light neutrino and gluino exchange mechanisms, and ii) 
\t)l\ 2 x 10 14 (thick dashed-dotted line) and |?7a'| 2 x 10 14 (thick dashed line) of the heavy LH 
Majorana neutrino and gluino exchange mechanisms, and using Ty 2 ( 76 Ge) = 2.23 x 10 25 
y and T^ 2 ( lm Xe) = 1.60 x 10 25 y. See text for details. 
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of these five nuclei, derived within the Self-consistent Renormalized Quasiparticle Ran- 
dom Phase Approximation (SRQRPA), were employed in our analysis. They correspond 
to two types of nucleon-nucleon potentials - Argonne ("Argonne NMEs") and CD-Bonn 
("CD-Bonn NMEs"), and two values of the axial coupling constant qa = 1.25; 1.00. Given 
the NMEs and the phase space factors of the decays, the half-life of a given nucleus de- 
pends on the parameters |?7 K | 2 of the mechanisms triggering the decay (eq. ([7])). 

In the present article we have considered in detail the cases of two non-interfering and 
two interfering mechanisms inducing the (/3/3) 0! ,-decay. If two non-interfering mechanisms 
A and B cause the decay, the parameters 1 77^4 1 2 and |?7b| 2 can be determined from data on 
the half-lives of two isotopes, T\ and T 2 as solutions of a system of two linear equations. 
If the half-life of one isotope is known, say Ti, the positivity condition which the solutions 
\t)a\ 2 and |//s| 2 must satisfy, |?7a| 2 > and \f]B\ 2 > 0, constrain the half-life of the second 
isotope T 2 (and the half-life of any other isotope for that matter) to lie in a specific 
interval [T]. If A and B are interfering mechanisms, 1 77^4 1 2 and 1 77^ | 2 and the interference 
term parameter, zab = 2 cos oabIvaVbI which involves the cosine of an unknown relative 
phase oab of r\A and T)b, can be uniquely determined, in principle, from data on the 
half-lives of three nuclei, T lj2 ,3. In this case, given the half-life of one isotope, say T 1; the 
"positivity conditions" \r]A\ 2 > 0, \t]b\ 2 > and —1 < cosoab < 1 constrain the half-life 
of a second isotope, say T 2 , to lie in a specific interval, and the half-life of a third one, 
T3, to lie in an interval which is determined by the value of T\ and the interval of allowed 
values of T 2 . 

For all possible pairs of non- interfering mechanisms we have considered (light, or 
heavy LH Majorana neutrino, and heavy RH Majorana neutrino exchanges; gluino, or 
squark- neutrino, and RH Majorana neutrino mechanisms), these "positivity condition" 
intervals of values of T 2 were shown in [T| to be essentially degenerate if T\ and T 2 
correspond to the half-lives of any pair of the four nuclei 76 Ge, 82 Se, 100 Mo and 130 Te. 
This is a consequence of the fact that for each of the five single mechanisms discussed, the 
NMEs for 76 Ge, 82 Se, 100 Mo and 130 Te differ relatively little [H[l3j: the relative difference 
between the NMEs of any two nuclei does not exceed 10%. One has similar degeneracy 
of "positivity condition" intervals T 2 and T 3 in the cases of two constructively interfering 
mechanisms (within the set considered). These degeneracies might irreparably plague 
the interpretation of the (/3/3)oi/-decay data if the process will be observed. 

The NMEs for 136 Xe, results of calculations of which using the SRQRPA method are 
presented in the present article, differ significantly from those of 76 Ge, 82 Se, 100 Mo and 
130 Te, being by a factor ~ (1.3 — 2.5) smaller. As we have shown in the present article, 
this allows to lift to a certain degree the indicated degeneracies and to draw conclusions 
about the pair of non-interfering (interfering) mechanisms possibly inducing the (/3/3)oi/- 
decay from data on the half-lives of 136 Xe and of at least one (two) more isotope(s) which 
can be, e.g., any of the four, 76 Ge, 82 Se, 100 Mo and 130 Te considered. 

We have analyzed also the possibility to discriminate between two pairs of non- 
interfering (or interfering) (/3/3)ot,-decay mechanisms when the pairs have one mechanism 
in common, i.e., between the mechanisms i) A + B and ii) C + B, using the half-lives of 
the same two isotopes. We have derived the general conditions under which it would be 
possible, in principle, to identify which pair of mechanisms is inducing the decay (if any). 
We have shown that the conditions of interest are fulfilled, e.g., for the following two pairs 
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of non-interfering mechanisms i) light neutrino exchange (A) and heavy RH Majorana 
neutrino exchange (B) and ii) gluino exchange (C) and heavy RH Majorana neutrino 
exchange (B), and for the following two pairs of interfering mechanisms i) light neutrino 
exchange (A) and heavy LH Majorana neutrino exchange (B) and ii) gluino exchange (C) 
and heavy LH Majorana neutrino exchange (B), if one uses the Argonne NMEs in the 
analysis. They are fulfilled for both the Argonne NMEs and CD-Bonn NMEs, e.g., for the 
following two pairs of interfering mechanisms i) light neutrino exchange (A) and gluino 
exchange (B), and ii) heavy LH Majorana neutrino exchange (C) and gluino exchange 
(B). 

We have also exploited the implications of the EXO lower bound on the half-life of 
136 Xe for the problem studied. We have shown, in particular, that for all four pairs 
of non-interfering mechanisms considered and the Argonne NMEs, the half-life of 76 Ge 
claimed in [21] is incompatible with the EXO lower bound on the half-life of 136 Xe |2U| . 
If we use the CD-Bonn NMEs instead, we find that the result half-life of 76 Ge claimed in 
|24] is compatible with the EXO lower bound on the half-life of 136 ATe for values of the 
corresponding LNV parameters lying in extremely narrow intervals. 

To summarize, the results obtained in the present article show that using the (/30)o v - 
decay half-lives of nuclei with largely different nuclear matrix elements would help re- 
solving the problem of identifying the mechanisms triggering the decay. 
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